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Revisiting Ulysses Observations of Interstellar Helium 

Brian E. Wood^, Hans-Reinhard Miiller^, Manfred Witte^ 

ABSTRACT 

We report the resnlts of a comprehensive reanalysis of Ulysses observations of 
interstellar He atoms flowing throngh the solar system, the goal being to reassess 
the interstellar He flow vector and to search for evidence of variability in this 
vector. We find no evidence that the He beam seen by Ulysses changes at all 
from 1994-2007. The direction of flow changes by no more than ~ 0.3° and the 
speed by no more than ~ 0.3 km s“^. A global fit to all acceptable He beam maps 
from 1994-2007 yields the following He flow parameters: Vjsm = 26.08 ± 0.21 
km s-\ A = 75.54 ± 0.19°, 13 = -5.44 ± 0.24°, and T = 7260 ± 270 K; where A 
and (3 are the ecliptic longitnde and latitnde direction in J2000 coordinates. The 
flow vector is consistent with the original analysis of the Ulysses team, bnt onr 
temperatnre is significantly higher. The higher temperatnre somewhat mitigates 
a discrepancy that exists in the He flow parameters measnred by Ulysses and 
the Interstellar Boundary Explorer^ bnt does not resolve it entirely. Using a 
novel techniqne to infer photoionization loss rates directly from Ulysses data, we 
estimate a density of tt-hb = 0.0196 ± 0.0033 cm“^ in the interstellar medinm. 

Subject headings: Snn: heliosphere — ISM: atoms 


1. INTRODUCTION 

The global strnctnre of the heliosphere is determined in large part by the local ISM 
flow vector in the rest frame of the Snn. Older determinations of this vector, particnlarly 
those based on measnrements of interstellar He flowing throngh the inner heliosphere by 
the GAS instrnment on Ulysses, have recently been challenged by new measnrements of the 
He flow by the Interstellar Boundary Explorer (IBEX). The canonical Ulysses analysis of 


^Naval Research Laboratory, Space Science Division, Washington, DC 20375, USA; 
brian.wood@nrl.navy.mil 

^Department of Physics and Astronomy, Dartmouth College, Hanover, NH 03755, USA 
^Max-Planck-Institute for Solar System Research, Katlenburg-Lindau, Germany 



2 


Wittd fl2004f) suggested the following He flow parameters: Vism = 26.3 ± 0.4 km s A = 


75.4 ±0.5°, /5 = —5.2 ±0.2°, and T = 6300 ±340 K; where A and (3 are the ecliptic longitude 

2000 coordinates. In contrast, IBEX d a ta recently yielded th e 


and latitnde direction in 
following flow properties flBzowski et al.l 120121: iMobius et al.l 120121: iMcComas et al.N2012l) : 
Vism = 23.2 ± 0.3 km s-^ A = 79.00 ± 0.47°, (3 = -4.98 ± 0.21°, and T = 6300 ± 390 K. 
Significant inconsistencies exist between the Ulysses and IBEX measnrements, particnlarly 
for Vism and A. 


The lower Vism valne of IBEX is of particnlar interest, as it may imply the nonexistence 
of a bow shock aronnd the heliosphere. The Vism = 23 — 27 km s“^ interstellar flow speed 
happens to imply a fast magnetosonic Mach nnmber of Mp ~ 1, making the existence or 
nonexistence of a bow shock aronnd the heliosphere very mnch an open qnestion. Uncer¬ 
tainties in the strength and orientation of the interstellar magnetic held, Bism, represent 
one obstacle in determining whether a bow shock exists. The higher Bism is (and the more 
perpendicnlar to the ISM how), the lower Mp should be, and the less likely there is to be a 
bow shock. Bnt nncertainties in Vism are also important. 


For many years the best assessm ents of ISM veloc i ty snggested Vism ~ 26 km s“^, not 
only the Ulysses /GAS measnrements flWitte et al.lll996l: IWittell2004l) bnt measnrements from 
ISM absorption lines as well (Vism = 25.7 ± 0.5 km s“^; Lallement & Bertin 1992). With 
these relatively high valnes, heliospheric modelers favored Mp > 1, implying the existence of 
a bow shock. However, not only has the new IBEX i neasnrement called this i nto qnestion, bnt 
so has one reanalysis of ISM absorption line data flRedheld fc Linskvl 1200811 . which snggests 
Vism = 23.84±0.90 km s“^. The lower Vism = 23 — 24 km s“^ valnes reported recently have 


been enongh for many to ar gue that Mp < 1 should be preferred flMcComas et ah 


2012 


Zieeer et al.ll2013fl . However, IScherer fc Fichtneii fl2014j) argue that including He^ density in 


the calcnlation of Alfven speeds instead of jnst assnming a pnre p roton plasni a wonl d still 
snggest Mp > 1 even if Vism ~ 23 km s“^. On the other hand, IZank et al.l (120131) note 
that charge exchange processes may tnrn the bow shock into more of a “bow wave” even if 
Mp > 1. 


In any case, the existence or nonexistence of a heliospheric bow shock is one issne driving 
interest in whether the He flow vector of Ulysses or that of IBEX i s to b e preferred. Assnming 
that the ISM clond aronnd the Snn is nonrigid, iGrv fc Jenkin.sl (120141) have reanalyzed ISM 
absorption line data and infer clond kinematics more c onsistent with the Uly s ses ve locity 


being valid close to the Snn. iBen-Jaffel et al.l (120131) . iLallement fc BertanxI (1201411 . and 


Vincent et al.l (1201411 have recently made additional argnments in favor of the older Ulysse s 


vector, implying that the IBEX measnrements mnst somehow be in error. iFrisch et ahl (l2013l) 
propose a very different solntion. They propose that the local ISM flow actnally varies with 
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time, and that these variations are responsible for the Ulysses/YBWK. discrepancy. However, 
observations of interstellar H from the Solar Wind Anisotropies instrument on the Solar and 
Heliospheric Observatory and the Space Telescope Imaging Spectrograph inst rument on the 


Hubb le Space Telescope do not provide support for a variable ISM flow vector fIVincent et ah 

2ni4h . 


All this attention on the He flow vector provides motivation for taking another long look 
at the old U lysses/GAS data, in order to see if a new, independent analysis confirms the 
Wittel (120041) results, and to see if there is any e vidence for He flow variability in the Ulysses 
data that would support the iFrisch et ahl (l2013h interpretat ion favoring a variable ISM flow 
vector. Further justification is that while the IWittel (120041) analysis considered data taken 
from 1990-2002, later data acquired up until the end of the Ulysses mission in 2007 were not 
considered. 


Katushkina et al.l (120141) have recently analyzed a couple of Ulysses He beam maps. 


including one from 2007, demonstrat ing th a t the new IBEX He vector cannot reproduce 
either map. They also note that the IWittel (120041) vector might in principle be able to fit 
the IBEX data i f the temperature is ahowed to be 9000 K, consistent with more extensive 


IBEX analyses ( Bzowski et al.l I2OI2I: iMobius et ahl I2OI2I: iMcComas et ahl 12012 ). We here 


present our own comprehensive reanalysis of the full Ulysses/GAS data set, providing a new 
assessment of the He flow vector usi ng analysis techn iques developed indep endently from 


both the original Ulysses / GAS team (jWitte et al.lll996l) and the IBEX team (iBzowski et ah 
2012 : Mobius et ah 2012 ). 


THE ULYSSES BATA 


Launched in 1990 October, the primary mission of the Ulysses spacecraft was to study 
various interplanetary constituents like the solar wind, magnetic helds, radio waves, dust, 
solar and galactic cosmic rays, etc, outsid e the ecliptic plane in an orbit passing over the solar 
poles (jWenzel et al.lll992l: lMarsdenll200l[) . This orbit was accomplished with a gravitational 
assist from Jupiter in 1992 February, leading to an orbit nearly perpendicular to the ecliptic 
plane, with an aphelion near Jupiter’s distance of 5 AU, and a perihelion near 1 AU. Ulysses 
continued to make observations from this orbit until the mission ended in late 2007. The 


Ulysses/GAS instrument, described in detail bv lWitte et al.l (119921) . was designed primarily 
to study interstellar neutral He atoms flowing through the inner heliosphere. Neutral helium 
provides the best diagnostic of the undisturbed ISM flow due to the high abundance of 
helium, and due to the fact that helium is not greatly affected by charge exchange during 
its journey through the outer heliosphere, in contrast to hydrogen and oxygen, for example. 
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The first He observations were made by Ulysses /GAS during its c ruise phase out to 
Jupiter (1990-1992). These data were hrst studied bv IWitte et ahl fjl993l ). However, in this 
paper we will focus on data taken during the operational phase of the Ulysses mission, with 
Ulysses in its polar orbit after the Ju piter encoun ter. These He data are generally of higher 
quality than the cruise phase data flWittd l2004l) . with better sampling of the He beam. 
Additional advantages of the operations phase data will become very clear below. 


Successful detection of interstellar He required that the inflow velocity be high enough 
for the He atoms to exceed the particle detection threshold of the GAS instrument. With the 
Ulysses orbital plane being basically perpendicular to the ISM flow vector, the detectability 
of He depended mostly on the spacecraft velocity. Only when Ulysses was in the near-Sun 
part of its orbit was the spacecraft traveling rapidly enough so that the incoming He had 
sufficient energy to exceed the detection threshold. Thus, the Ulysses/ GAS He observations 
were conhned to the three fast latitude scans made by Ulysses, in the periods 1994 -1996, 
2000- 2002, and 2006-2007. These are the three epochs that we will focus on here. 

(|2004j) considered the hrst two epochs of observations, but not the last epoch. 


Witte 


The GAS device works essentially like a pin-hole camera with a single-element detector, 
so observing the He beam requires making many different exposures in different directions to 
fully map out the count rate distribution. The Ulysses spin axis, which was always pointed 
towards Earth, dehned how this was done. Ulysses had a rotation period of 12 s. During 
this period, the GAS instrument would observe at a particular elevation angle, measured 
relative to the spin axis, mapping out a ring on the celestial sphere. Detected counts were 
accumulated in 32 bins, which could uniformly cover the full 360° rotation; or more commonly 
cover smaller arcs of 90°, 45°, or 22.5°. After a period of typically 68 min, the elevation angle 
scanned by GAS would be changed by a step of 1°, 2°, 4°, or 8°; allowing another arc to be 
observed. Over the course of 2 — 3 days, typically, a full scan across the He beam would be 
completed, and a new map begun. 


The GAS instrument actually utilized two different channels with independent detectors, 
one with a wide held of view (WFOV) and one with a narrow held of view (NFOV), and 
both channels were exposed simultaneously. Over the course of the three fast latitude scans 
mentioned above, Ulysses obtained ~ 400 WFOV maps and ~ 400 NFOV maps. We will 
be working almost exclusively with the WFOV data here, which have higher signal-to-noise 
(S/N) than the NFOV maps. 


Figure 1 shows a typical WFOV map, from 2001 January 24. The GAS instrument 
suhers from some sensitivity to UV emission even when in its particle detection mode 
fjWitte et al.lll996l ). so some hot stars are visible in the map, in addition to the broad He 
beam. The stars can be inconvenient if blended with the He beam, but are actually useful 
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as calibrators of pointing accnracy otherwise. In this article we will in fact be nsing a re¬ 
cent recalibration of the Ulysses j GAS data set nia ^ng extensive n se of the stars as pointing 
calibrators, which is described in more detail bv iBzowski et ahl (120141) . In an end-to-end 
test, the calibration of the direction determination was adjusted so that the positions of 
stars observed by the instrument coincided with their astronomical position. Based on many 
more star observations from the whole mission a minor further adjustment of a few tenth of 
a degree was required and subsequently all He data were reprocessed in 2013 with this new 
calibration. This recalibration provides yet another motivation for reanalyzing these data, 
though the change in the data is modest. Note that although the archival Ulysses data, and 
its original analyses reported in the Witte et al. series of papers, are in B1950 coordinates, 
the numbers and hgures we provide in this article will be in the more modern J2000 system. 

A hnal distinction to be made is that we will only be considering observations of the 
so-called “direct beam” formed by He atoms that flow directly to the detector (with modest 
gravitational deflection), as opposed to the “indirect beam” created by particles that go 
around the Sun and are gravitationally redirected back towards the detector, arriving from 
a very different direction than the direct beam atoms. The indirect beam is much fainter 
than the direct beam due to massive photoionization and electron-impact ionization losses 
that occur when these particles pass close to the Sun. Ulysses has nevertheless detected the 
indirect beam in observations from 1995 fjWittel 120041) . but we will not be considering those 
data here. 


3. EMPIRICAL ANALYSIS 

Rather than immediately try to extract a He flow vector from the data, which requires a 
somewhat complex forward modeling analysis (see Sections 4 and 5), we hnd it worthwhile to 
hrst make some purely empirical investigations of the beam. The primary motivation for this 
is to look for any evidence that the beam is changing at all during the 1994-2007 time period. 
We quantify the properties of the beam using 2-D Gaussian hts to the beam, as mapped 
into ecliptic coordinates (as in Figure 1). Note that there is no physical reason the beam 
should be precisely Gaussian, but in practice the beams observed by Ulysses are reasonably 
Gaussian in shape, so this is a practical way to quantify them. An example Gaussian £t is 
shown in Figure 2. There are six free parameters of such a £t: the central longitude and 
latitude, the latitudinal and longitudinal widths [specihcally the full-width-at-half-maximum 
(FWHM)], the count rate amplitude at beam center, and a background count rate; where 
we assume a flat background underneath the He beam. 

In the following discussion, we make use of an orbital phase for Ulysses, cj), based on the 
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angle defined by the positions of the spacecraft and the Snn, and the location of the ecliptic 
plane crossing near aphelion. In this system, the ecliptic plane crossing near aphelion at 
ecliptic coordinates (A,/3) = (157.4°,0°) is at 0 = 0.0, and that near perihelion at (337.4°,0°) 
is at 0 = 0.5. The 0 = 0.25 and 0 = 0.75 orbital phases correspond to the near sonth pole 
and near north pole locations (relative to the Snn), so Ulysses is below the ecliptic at 0 < 0.5 
and above it at 0 > 0.5. (Ulysses travels from sonth to north dnring its fast latitnde scans.) 

The Ganssian fits are performed using a semi-automated procedure, but all the fits are 
visually verified. Many beam maps end up being excluded for various reasons. At 0 < 0.2 
and 0 > 0.9 the data are simply deemed too noisy. Observations with 0 < 0.26 are excluded 
due to contamination from a star (specihcally the rightmost star seen in Figure 1). Some 
maps are excluded for having a background that is not flat, possibly due to solar activity at 
that time. The most common reason for exclusion is simply that some scans did not fully 
cover the He beam. We only want to consider beam maps where the scan fully encompasses 
the He beam. Ultimately, fits are accepted for 234 NFOV maps and 238 WFOV maps, 
including the map in Figure 1. As mentioned above, we will be focusing on the higher S/N 
WFOV data, but we always do verify that results are not different when considering the 
NFOV data instead. 


The Gaussian £t parameters of the WFOV maps are plotted in Figures 3 and 4. Fig¬ 
ure 3a shows that the He beam traces out a horseshoe pattern in the sky during a fast latitude 
scan. Performing polynomial fits to longitude and latitude as a function of orbital phase al¬ 
lows us to define the average beam location illustrated in Figure 3a. We will emphasize 
below how the shape of this horseshoe is a powerful diagnostic for the He flow parameters. 
Figure 3a demonstrates that basically the same horseshoe pattern is made during all three 
fast latitude scans, showing little indication for variation in beam location. Figures 3b-d 
provide a better indication for the deviations of beam position about the average. The scat¬ 
ter of the data points is indicative of the measurement uncertainties, but some systematic 
deviations are discernible despite the scatter. For example, the beam longitude appears to 
be lower during epoch 3 by a few tenths of a degree. The beam latitude is also a few tenths 
of a degree lower during epoch 3, but only for 0 = 0.5 — 0.7. The latitudes observed during 
epochs 1 and 2 are discrepant for 0 = 0.75 — 0.85. 


Of particular interest are longitude variations, as it is in longitude that the IBEX and 
Ulysses flow vectors are in disagreement, and it is in longitude that iFrisch et al.l (120131) 
claim evidence for a systematic increase with time. There is, however, no evidence for any 
increase in longitude during the Ulysses era, with the epoch 3 data suggesting if anything 
a slight decrease, as noted above. The average discrepancy in Figure 3d is 0.28°. This can 
be interpreted as an estimate of the pointing uncertainties inherent in the Ulysses data. 
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consistent with expe ctations based on the pointing calibration done nsing stars fjWittell2004J: 


Bzowski et al.l 120141) . 


In any case, the longitnde variations in Figure 3b are much smaller than the 3.6° dis¬ 
crepancy in longitude between the IBEX and Ulysses He flow vectors. We conclude that 
there is no indication from the Ulysses data to support the notion that He flow variability is 
responsible for the IBEX/ Ulysses discrepancy. Furthermore, given that we now know that 
the 2007 Ulysses data are not more consistent with the IBEX results than the older data, 
there is only a 2 year time gap between the Ulysses and IBEX eras for the He flow to have 
changed so signihcantly. 


Figures 4c-d also demonstrate that there is no evidence for any signihcant change in 
the He beam shape from 1994-2007. The latitudinal and longitudinal beam widths change 
with orbital phase, but the same sort of variation is observed during all three fast latitude 
scans. Figures 4a-b show that the three epochs do, however, exhibit noticeable differences 
in background and beam amplitude. The Ulysses/GAS background is a complicated mix of 
solar and cosmic sources, and we will not discuss it much here, other than to note that both 
orbital phase and epochal variations are indicated by Figure 4b. 


As for the beam count rate amplitude in Figure 4a, it naturally peaks near </> = 0.5 
during each fast latitude scan, as that is when the spacecraft is moving fastest. The different 
count rate levels observed during the three epochs are most naturally interpreted as being 
indicative of variations in the solar EUV photoionization rate. It makes perfect sense for 
epoch 2 (2000-2002) to have the lowest He fluxes, as this is a time near solar maximum, 
with high photoionization rates and therefore substantial losses of neutral He due to pho¬ 
toionization. In contrast, epochs 1 and 3 are both near solar minimum, thereby minimizing 
photoionization losses. We will return to the related issues of photoionization and He density, 
riHe, in Section 5. 


4. MEASURING THE ULYSSES He FLOW VECTOR 

The simple quantihcation of the He beam properties provided by the Gaussian htting 
in the previous section is not helpful for actually inferring a He flow vector from He beam 
maps, which instead requires a forward modeling approach. In this procedure, a Maxwellian 
He velocity distribution is assumed at inhnity, characterized by hve parameters ( riHe , U ism, 
A, /3, T), and then this Maxwellian is propagated to the position of the spacecraft. The 
distribution will be primarily affected by the Sun’s gravitational influence. Loss processes 
can also be accounted for in this forward modeling, the most important being photoionization 








by solar EUV photons. We will initially ignore photoionization in our forward modeling, but 
we will consider it later to see if it changes our inferred He flow parameters. 


The f orward modeling method ut ilized in the original Ulysses /GAS team analyses is de¬ 
scribed by Banaszkiewicz et ah 1 199^. Our forward modeling routine, w' hich is described in 
detail elsewhere f Muller fc Cohen 2012 : Muller 2012 : Muller et ah 20131) . involves a rather 
different approach to the problem than the usual method involving numerical integration 
from infinity. In our computations, quantities that are conserved along a trajectory, such 
as total energy, angular momentum, eccentricity, and direction of perihelion, are used to 
determine in one algebraic step the neutral He velocity at a desired point in space when 
its velocity at inhnity is given. Conversely, when both location and velocity of a He atom 
are given, the conservation equations dictate a unique corresponding velocity at inhnity. 
The entire phase space distribution, F(v), can be calculated at a desired point in the inner 
heliosphere by decomposing the problem into a family of trajectories with different veloci¬ 
ties (different conserved quantities) that all converge at the chosen point. By assuming a 
thermal Maxwellian as the original distribution at inhnity, each of those trajectories can 
then be assigned its corresponding phase space density, yielding F(v) at the point of interest 
f Muller fc Cohen 2012 : Muller 2012 : Muller et ah 2013 1. 


Computing count rates in a given direction involves hrst transforming F(v) into the 
spacecraft rest frame and converting into a coordinate system dehned around the view direc¬ 
tion of interest, resulting in F{Vs,6,'ijj), where Vg is the speed relative to the spacecraft, 6 is 
the poloidal angle with respect to the view direction, and i/ is the azimuthal angle. With a 
volume element = Vg sin 6 dVg dip dO" , the predicted count rate can then be expressed 
as 


C = A Vs De{Vs) G{e) F{Vs,9,iP) d\, 


( 1 ) 


where A = 0.0908 cm^ is the ehective area, D^iVs) is the detection efficiency, and G{6) 
is the detector geometry function. F igure 5 shows our di g itized versions of G{9) and He, 
approximated from Figures 1 and 3 of lBanaszkiewicz et al.l fll996l ). with the latter expressed 
as a function of incoming particle energy. Eg = ^mVg, with m being the mass of the He 
atom. For a given elevation angle, particles in the azimuthal bins described in Section 2 are 
observed while the spacecraft is spinning, so there is “scan smoothing” along the azimuthal 
direction. We account for this by averaging G values computed at pixel center and at ±Aa/2, 
where Aa is the azimuthal angular bin size. 


By computing G for all n observed directions, we compute a synthetic He beam map, 
Gi (for i=l,n), for comparison with the real one. This synthetic beam is multiplied by 
a normalization factor, N, in order to best match the count rates of the observed beam. 
This amounts to a revision of the assumed riHe value to njjg = N x riHe- In making this 
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normalization we are in effect replacing riHe with N as the true density-related free parameter 
of our £t. We also have to account for the background count rate, B. As we did in the 
Gaussian fitting, we simply assume a flat background. Thus, the new predicted count rates 
are 


C' = NxCi + B. 


( 2 ) 


We perform a minimization (e.g., iBevington fc Robinsonlll992h to determine the N and 
B values that lead to the best fit to the data. The value of this best fit represents our 
measure of merit for this particular synthetic He beam map, which can be expressed as 


X 


2 



( 3 ) 


where C[ = G'[A^, 5, A,/?, T)] are the predicted count rates. Si are the observed 

count rates, cTj their Poissonian uncertainties, and n the number of points considered in the 
He beam map. We then have to determine what He flow parameters {Vjsmi A, 13, and T) 
truly minimize representing the best £t to the data. 


The standard approach for determining a set of best-fit parameters is to start with a 
guess for the parameters, compute x^) and then use a x^ minimization routine to refine 
the parameters to reduce x^ until ultimately the true x^ minimum is reached. If u is the 
number of degrees of freedom of the fit (the number of data points minus the number of free 
parameters), then the reduced chi-squared is defined as xl = whi c h sho uld be ~ 1 

for a good fit. We use the Marquardt method described by IPress et al.l fjl989l) as our x^ 
minimization routine. 


Figure 6a sho ws a WFOV map from 1995 August 20, and Figure 6b-c shows a fit to 
this map, with the IWittel (120041 ) vector used as our initial guess for the ht parameters. It 


is a simple m atter for the particle tracking code mentioned above flMiiller fc CohenI 12012 


Mlillerl 120121) to compute the expected beam location based on the IWittel (120041) vector, if 


we assume T = 0 K (meaning all He atoms have the exact same velocity and follow the 
exact same trajectory). We use this prediction to define the map region to be fitted. In 
particular, we only consider data points within 25° of this location (white dots in the figure), 
which encompasses enough of the region beyond the He beam to estimate the background. 
Figure 6b shows the synthetic beam map of the best fit, and Figure 6c explicitly plots 
observed versus predicted count rates, with a data point plotted for each of the white dots 
in Fig u re 6b . The He flow vector parameters listed in Figure 6c are gratifyingly close to the 


Wittel (120041) vector, and the xl = 1-06934 value indicates an excellent quality fit. 


Unfortunately, the best-fit parameters that we infer are actually highly dependent on 
our initial guesses for the parameters. This is shown explicitly in a second £t to the 1995 An- 
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gust 20 data in Figure 6d-e. With an initial guess very different from the IWittel (12004) vector 
used in Figure 6b-c, we end up with a fit to the data with very different parameters. The 
parameters listed in Figure 6e are very far from all modern estimates of the local ISM flow 
vector, so we know they are wrong. But these parameters lead to a fit that is just as good as 
the more plausible fit in Figure 6b, with & xt — 1-05731 value that is actually slightly lower 


even. 


This exercise emphasizes just how difficult it is to find a unique He flow vector from 
a single He beam map. A single map does not sufficiently constrain the problem. There 
are significant degeneracies among the fit parameters, leading to long diagonal troughs in 
space where it is very hard to find the true minimu m along the trough. This is true 
for the analysis of IBEX data as well (see Figure 22 in iBzowski et al.ll2012l) . The most 
basic degeneracy is between flow direction and flow speed, which can be explained by the 
following question: If a beam is seen from a given direction, is that due to a fast beam 
coming from roughly that direction, or due to a slower beam that has been gravitationally 
deflected into that direction? Since the He flow vector happens to he near the ecliptic plane, 
this degeneracy basically ends up mostly between Vjsm and A (rather than /3). 

We argue here that these degeneracies can be effectively broken if Ulysses data from 
different parts of the Ulysses orbit are fitted together, as opposed to fitting one beam map 
at a time, as in Figure 6. The utility of this approach is illustrated by Figure 7. We had 
previously emphasized how the He beam traces out a horseshoe pattern in the sky during a 
fast latitude sca n (see Figure 3a). Figure 7 shows the horseshoe pattern expected based on 
the IWittel (12004) vector (assuming T = 0 K), and also shows what happens to the horseshoe 
when /?, A, and Vjsm are changed. When (3 and A are changed the horseshoe shifts vertically 
and horizontally, respectively. In contrast, Vjsm affects the size of the horseshoe, with higher 
ViSM yielding a smaller horseshoe. The crucial point is that Vjsm affects the horseshoe shape 
very differently from f3 and A, thereby breaking the direction/velocity degeneracy noted 
above. But this diagnostic power can only be properly considered if He beam maps from all 
parts of the horseshoe are considered simultaneously in fitting the data; hence the global fit 
approach. 


The specific data used in this global fit approach are the 238 WFOV maps used in 
the Gaussian fitting (see Section 3). The global fit is similar to the individual map fits in 
Figure 6, but we are now fitting many more data points, naturally; 81,405 count rates within 
the 238 maps, to be precise. In addition to the four He flow parameters of interest (Wsm, 
A, /?, and T), there are also N and B parameters for each individual map, for a total of 480 
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free parameters. Summing over the 238 separate maps, equation (3) becomes 


X 


238 n 

j=i 1=1 


Sj,i 




a. 




(4) 


where C' j = C'JiV,-, A,/9, T)]. Figure 8 shows the 238 Nj and Bj values that 

result from our best £t. The Bj values are essentially identical to the backgrounds shown 
in Figure 4b, measured in the Gaussian fit analysis. The Nj factors are directly related 
to inferred He densities (see above), so that is how they are shown in the hgure. The clear 
phase dependence is an effect of photoionization losses, which are more severe at intermediate 
phases when the spacecraft is closer to the Sun. We will discuss this issue in detail in Section 5 
when we actually seek to correct for photoionization. 


We need to not only find the best fit but also to estimate uncertainties in the fit param¬ 
eters, so the full analysis actually involves a series of fits with one of the four flow parameters 
held constant, as we scan through the parameters to see how varies with Vism, A, /?, and 
T. Results are shown in Figure 9. There is a very clear minimum, Xmim each panel. We 
define = X^~Xminy each panel of Figure 9 shows the variation of Ay^ across the Xmin 
region. Fifth order polynomials are htted to the data points to int erpolate between them. 
The A y^ v alues are used to d efine the error bounds around Xmin- iBevington fc Robinson 


(I 1992 I ) and IPress et al.l (1l989l ) both provide useful discussions about how best to do this. 
For the number of free parameters of our fi t (480), the 3a confidence con tour corresponds 
to Ay^ = 578.3, based on relation 26.4.14 of lAbramowitz fc Steeunl (1196511 . This is the Ay^ 
level used to define the uncertainty ranges shown in Figure 9. 


Our final derived He flow parameters are: Vjsm = 26.08 ± 0.21 km s“^, A = 75.54 ± 
0.19 deg, (3 = —5.44 ± 0.24 deg, and T = 7260 ± 270 K; which are listed in the first line of 
Table 1. Our best global fit has xl = 1.524, which is somewhat higher than it should be 
ideally, indicating a modest level of systematic discrepancy from the data. It is not practical 
to show the actual fit to the data for all 238 beam maps, but for illustrative purposes 
Figure 10 shows the fit for two particular maps, one from 1994 November 9 and another 
from 2007 November 26. The former exhibits no clear systematic deviation between the fit 
and the data, but the latter shows some systematic deviation. 


Calculations involving photoionization will be discussed in more detail in Section 5, but 
at this point we mention that in order to see whether photoionization could be affecting our 
He vector and T measurements, we tried redoing the global fit assuming a He photoionization 
rate at 1 AU of /3inr7 = 1.5 x 10~’^ s~^, which i s near the high end of those typically estimated 
( Bzowski et al. 2012 : Katnshkina et al.ll2ni4li . The resulting parameters and y^ are: Vjsm = 
26.11±0.21 km s-S A = 75.47±0.18 deg, f3 = -5.49±0.23 deg, T = 7330±280 K, xl = 1-528. 
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These values are very similar to those found before, so considering photoionization does not 
change the derived vector or t emperature signihcantly, and does not affect the quality of £t, 
consistent with the findings of iKatushkina et ahl (120141) . 


Figure 9 and Table 1 demonstrate that our results ar e gener a lly in very good agreement 
with the previous canonical Ulysses measurements from IWittel (120041) . with the exception 
of T. Our t emperature is significantly higher than the T = 6300 ± 340 K measurement of 


Wittel (120041) ■ As discussed in Section 2, we are working with a new reduction of the Ulysses 


data, but we have performed fits to He beam maps from the older data reduction, and we do 
not hnd any signihcant effect on our derived temperature. We initially suspected that our 
higher T measurement may be an artifact of the global fitting approach. To first order, T is 
determined mostly by the size of the observed He beams, while the other three parameters 
determine its direction. If the global £t has difficulty simultaneously reproducing the beam 
center for all 238 beam maps, it will try to compromise by broadening the beam, increasing 
the inferred T. T his effect will not be an issue with single-map fits, and so would not affect 
the IWittel (1200411 analysis approach. 


We perform single-map fits to our 238 beam maps, where we fix Vjsm = 26.08 km s“^ 
in order to resolve the parameter degeneracy problem with single-map fits described above. 
This provides 238 measurements of A, (3, and T. The average and standard deviation of these 
values are reported in the second line of Table 1. The A and (3 measurements are nearly 
identical to the global fit values. The T = 7090 ±370 K measurement is somewhat lower than 
the T = 7260 ± 270 K global fit value, as expected, but not by that much. Experimenting 
with other velocity values within our quoted Vjsm = 26.08 ± 0.21 km s“^ range does not 
change things significantly. We conclude that the beam-broadening effect of the global fit is 
not the dominant cause of our higher T. 


Differences in background treatment are another possible cause of the T discrepancy. 
The original Ulysses data analyses exclu ded pixels with count rates too close to the back¬ 
ground value (jBanaszkiewicz et al.Nl996l) . Typically only the upper ~ 80% of the prohle 
was considered, whereas our approach is to consider all pixels within 25° of the beam, with 
the background determination actually being part of the fit rather than assuming a prede¬ 
termined background. We have considered the possibility that focusing more on the central 
part of the beam could lead to the perception of a narrower beam and therefore a lower 
T. However, using single-map fits we have investigated the effects of focusing on the cen¬ 
tral part of the beam and we do not find that such fits consistently lead to significantly 
lower temperatures. Thus, we ultimately dismiss this possible explanation for our high T 
measurement. 


We note that our higher temperature is in better agreement with the T = 7500 ± 
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1300 K Local Interstella r Cloud (LIC) temperature estimated from ISM absorption lines 
flRedfield fc Linsk'vll2008l). It is also coii s istent with another contemporaneous reassessment 


of the Ulysses He data bv iBzowski et al.l (120141). largely follo wing the same analysis approach 
used previously to study IBEX data flBzowski et alJl2012l) . Their derived best-£t He flow 


parameters (Vism = 26.0 km s A = 75.3°, {3 = —6°, T = 7500 K) are in reasonably good 
agreement with our results, including the higher temperature. 


The only difference is that iBzowski et ahl (1201411 claim much larger uncertainties. Specif¬ 
ically, the parameters are quoted as having uncertainty ranges of 74.2° < A < 76.5°, 
-7° < /3 < -5°, 24.5 < Vism < 27.0 km s-\ and 5500 < T < 9000 K. The large er¬ 
ror bars are based on independent analyses of data taken from different parts of the Ulysses 
orbit, but we consider these error bars to be overly conservative, as this analysis approach 
greatly weakens the diagnostic power of considering data from all parts of the Ulysses orbit 
simultaneously, which is crucial to the efficacy of the global £t approach. Figure 7 indicates 
how A, 13, and Vjsm change the shape of the horseshoe-shaped track for the He beam in 
different ways, and it is this effect that allows the global £t approach to break the degen¬ 
eracies that plague single-map analyses (see Figure 6). Considering observations from only 
one part of the Ulysses orbit will in effect reintroduce the degeneracy problem , naturally 
leading to the inference of the large uncertainties quot ed bviBzowski et al.l (20141) . To put it 
another way, we see the large uncertainties quoted by Bzowski et al. ~ (l2014 ) as being a fur¬ 
ther demonstration of the importance of considering data from all parts of the Ulysses orbit 
together, rather than being an accurate quantification of systematic uncertainties inherent 
in these data. 


The cause of the higher T inferred both by ourselves and iBzowski et ahl (120141) remains 
to be explained. Having excluded the po ssibilit i es me ntioned above, we conclude that the 


most likely cause of this discrepancy with IWitte (l2004f) lies in t he Ba yesian statistics of the 
older analyses. While our analysis and that of IBzowski et al.l (120141) rely solely on the 
statistic to quantify how well a set of parameters is htting the data, the older analyses instead 
use a Bayesian merit function that considers not only consistency with the data, but also 
consi stency with prior assumpti ons about what the best-fit parameters should be (see eqn. 10 
from iBanaszkiewicz et al.lll996[) . The advantage of this Bayesian approach is that this is one 
way to break the degeneracies in single-map fits, a nd derive rela tively self-consistent single¬ 
map measurements like those shown in Figure 2 of IWittd (120041) . The disadvantage is that 
the analysis is potentially biased by the a priori assumption of expected best-£t parameters, 
as well as the question of how to weight the relative importance of consistency with those 
parameters and consistency with the data. The clearest distinction between our analysis and 
the older ones lies in this difference in how the best £t is dehned, so we propose that the 
Bayesian approach of the older work drove the analysis towards a lower temperature, while 
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the pure x^-minimization analysis used here and by iBzowski et ahl (120141) leads to a higher 
temperature. 


The higher temperature suggested by both our analysis and that of iBzowski et ahl (120141) 


could have relevance for the discrepancy in the IBEX and Ulysses He flow vectors. The IBEX 
analysis suffers from the same kind of parameter degeneracies di scussed above for U lysses. 
Although the IBEX data suggest a minimum for T 6200 K (IBzowski et ahl 120121 ). it is 


not clear that significantly higher temperatures are excluded by the data. Furthermore, the 
parameter dependencies for IBEX are such that assuming a higher temperature would push 
the IBEX-derived A and Vjsm values closer to those of Ulysses. Specihcally, assuming our 


and Vtsm = 24.6 ± 0.4 km s 


-1 


Bzowski et ah 

2012 ; 

Mobius et ah 

2012 : 

McComas et ah 


20121) . These values are still n o t qui te consistent with Ulysses, but they are signihcantly 
closer than the iMcComas et ahl (120121 ) values quoted in Section 1. 


A likely important source of systematic uncertainty in the Ulysses data analysis that 
we want to emphasize is uncertainty in the detection efficiency (He) curve in Figure 5b. It 
is important to note that He atoms are typically detected in the = 30 — 60 eV range. 
This is very much in the “knee” of the curve, where D^, is signihcantly energy dependent. 
Combine this energy dependence with the fact that there is generally a signihcant gradient 
of average particle energy across a He beam, and the result is a situation where particle 
detection efficiency is higher on one side of the beam than the other, leading to an observed 
beam that is shifted from where one might naively expect to see it. This is clearly the case 
in Figure 6 for example, where the observed beam is shifted to the l eft fro m its expected 
location based on a zero-temperature calculation assuming the IWittel (120041) vector. 


This effect is shown in more detail in Figure 11, which compares the observed trace of 
the He beam center from Figure 3a with a trace computed from our best-fit He ffow vector, 
assuming T = 0 K. The average discrepancy between these two horseshoe-shaped tracks is 
1.40°, which is at least mostly due to the effects of described above. The magnitude 
of this effect will be quite sensitive to the slope of the He curve around Eg = 45 eV. We 
therefore suspect that uncertainties in this slope are an important source of systematic error 
in the Ulysses analysis. 


Could this effect in principle resolve the discrepancy in the IBEX and Ulysses He flow 
vectors? To investigate this possibility. Figure 11 shows the location of t he ho rseshoe track 
that Ulysses should have observed if the IBEX vector of lMcComas et ahl (120121) was correct, 
approximated once again with the T = 0 K assumption. In most places along the horseshoe, 
the He slope effect is not shifting the track towards the IBEX track, so assuming Hg curves 
with different slopes will unfortunately not resolve the IBEX/ Ulysses discrepancy. The IBEX 
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track in Figure 11 is an average of 3.73° from the observed track, illustrating the magnitude 
of the disagreement in angular terms. Based on Figure 7, the lower Vjsm and higher A values 
measured by IBEX (see Section 1) would be expected to produce a horseshoe track that is 
broader and shifted to the right from that observed by Ulysses, and Figure 11 shows that 
explicitly. 


Finally, we revisit the issue of He flow variability by performing global hts to the Ulysses 
data separately for the three fast latitude scans, leading to independent measurements for 
1994-1996, 2000-2002, and 2006-2007. Results are listed in Table 1 and displayed in Fig¬ 
ure 12. Consistent with the Endings of the empirical analysis of Section 3, we see no evidence 
for any He flow parameter variability, with the He flow direction varying by no more than 
~ 0.3° from 1994-2007, and the velocity varying by no more than ~ 0.3 km s“^. With 
IBEX data acquisition beginning in early 2009, there is less than a two year time gap for the 
3.5° jump in flow longi tude to occur, if th e variable flow interpretation of the IBEX/ Ulysses 
discrepancy is correct (iFrisch et al.l 1201311 . 


DENSITY AND PHOTOIONIZATION 


The previous section described how we measured four of the five parameters describing 
the assumed Maxwellian He distribution in the ISM. The remaining parameter is the density, 
riHe- Unlike the other four parameters, the density is significantly affected by loss processes 
that occur during the He atoms’ journey from the ISM to Ulysses. Thus, measuring riHe 
is somewhat more complicated. The loss rate, (3ion(j), is assumed to vary as r~^, and is 
generally quantified by quoting the loss rate at 1 AU. For He, the most dominant loss 
process is photoioniza tion by solar EUV photons, thou gh electron impact ionization may 
also be non-negligible fjWittell2004 iBzowski et al.ll2013l) . For an assumed (3ion at 1 AU, our 
particle tracking code can account for losses along the track, assuming the dependence. 
At the ISM flow speed it takes about 500 days to travel from 10 AU inward to the observer 
(see, e.g.. Figure 3 in Witte 2004), so the effective ionization rate should be an average of 
the rates encountered over the past year or so, weighted as r“^. 

The usual procedure for dealing with pho toionization is to rely on estimates of R 


from various moni t ors of solar EUV emission dMclV 


Woods et al. 

2005; 

Lean et al. 

2011 : 


Bzowski et ah 


Vlullin et al. 

2004: 

2013 

). However, 


significant uncertainties, as many of the EUV irradiance proxies are less than ideal for 
our purposes, and some of them are inconsistent. The Solar EUV Experiment (SEE) on 
the Thermospheric Ionospheric Mesospheric Energy and Dynamics (TIMED) spacecraft has 
significantly improved matters, but it was launched in 2002, and so does not cover most of the 
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Ulysses time period. iLean et al.l (120111) quote uncertainties in EUV irradiance of 10%-20% 
in the SEE era, but up to a factor of 2 uncertainties before it. The SEE EUV measurements 
have been correlated with other solar activity indicators that have been monitored for many 
decades (e.g., Mg II, E 10 . 7 ), allowing those indicators to be utilized as EUV proxies. This 
has improved matters even for the older time period. 


Nevertheless, we here try a different technique for dealing with photoionization that uses 
only Ulysses constraints. Focusing on the second fast latitude scan (2000-2002), Figure 13a 
shows densities measured as a function of orbital phase, computed assuming different pho¬ 
toionization rates. When an overly low rate is assumed (/djon = 0 s“^ in the hgure) the result 
is a U-shaped curve with low densities at intermediate phases (0 = 0.4 — 0.6). Although 
densities at all phases will be underestimates, the effect is more pronounced at these interme¬ 
diate phases, since that is when Ulysses is closest to the Sun and therefore observes a He flow 
with the highest losses. Conversely, when an overly high rate is assumed (/djon = 3.0 x 10“^ 
s“^ in the figure) the pattern reverses and overly high densities are observed at intermediate 
phases, indicative of an overcompensation in the loss correction. 


We define the most likely photoionization loss rate as that which leads to the flattest 
density-vs-phase curve. Specihcally, we compute the standard deviation of the density mea¬ 
surements divided by the average density, o'n{,f^ion)/n{f^ion), and we find the (3ion value that 
minimizes this quantity. At this minimum, the average density n is then our riHe mea¬ 
surement, and the minimum deviation is (Tn,min/n-Re- We estimate la uncertainties for our 
loss rate and density measurements using the standard deviation of the mean at this mini¬ 
mum, amean = o'n,minlV^m, wliere here Nm is the number of measurements. The range of 
acceptable loss rates are then defined to be where 






tonj ^ ^n,min “ 1 “ 


^He 


( 5 ) 


For the epoch 2 data in Figure 13a, the minimum density variance is at f3ion = 1-79 x 10“^ 
s“^, and the resulting densities are shown in the hgure. Assuming this rate, the average 
He density is then computed to be nne = 0.0197 cm“^. Figure 13b displays the density 
measurements for all three epochs. The inferred densities and photoionization rates are 
listed in Table 2, along with the la uncertainties computed as described above. For epoch 3, 
the analysis is less constrained due to fewer data points and intrinsically greater scatter (see 
Figure 13b), so we essentially hnd only an upper limit to (3ion- 


The densities independently measured from the three epochs are not identical. Assum¬ 
ing that the ISM He density does not really vary, the degree of inconsistency is presumably 
an indication of the uncertainties in the loss correction, although detector sensitivity uncer¬ 
tainties are also a potential factor. The stability of the detector sensitivity was tested twice 
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during the course of the mission. New layers of LiF were evaporated onto the surface from 
the GAS instrument’s LiF supply in 1995 and 2001. The He signal was measured immedi¬ 
ately before and after each deposition. It showed that the detection efficiency did not change 
within the available statistical accuracy of 20%. Unfortunately, a third test planned for 2006 
had to be cancelled, as the power necessary to heat the LiF supply may have exceeded the 
power available to the spacecraft at the time, which had decreased during the course of the 
mission. Thus, it cannot be proven that the detection efficency did not degrade during the 
third and hnal epoch, altough it seems unlikely in view of past experience. In any case, the 
average and standard deviation of the three measurements in Table 2 provides our hnal best 
estimate for the ISM He density: nne = 0.0196 ± 0.0 033 cm~^. T his is somewhat higher 
than the nne = 0.015 ± 0.003 cm“^ measurement from IWittel (120041) . due at least in part to 


higher assumed photoionization rates, but the error bars overlap. 


T able 2 compares our inferred loss rates with average ones estimated by iBzowski et ah 


(120131. 120141) for the three epochs based on solar wind and EUV constraints, taking into 
consideration both photoionization and electron impact ionization. There is actually good 
agreement for epochs 2 and 3, but not for epoch 1, where our jSion = (1.36 ±0.27) x 10“^ s“^ 
is signihcantly higher than the (3ion = 0.9 x 10“^ s“^ rate expected for this period. Given that 
our estimated density is also rather high, it seems likely that our analysis approach has led 
to an overestimate of /Sjon- Our technique of minimizing phase-dependent density variation 
does have the drawback of not accounting for any change in photoionization rate that might 
occur in the 20 month peri od covering the (j) = 0 . 25 —C ).90 fast-latitude scan range. However, 
the loss rates estimated by iBzowski et ahl (l2013l. 120141) do not actually show much variation 
during epoch 1, so unidentihed issues may also be present. 


SUMMARY 


Motivated in part by the discrepancy in the He flow vectors measured by Ulysses and 
IBEX, we have performed an independent analysis of the Ulysses/GAS data set from 1994- 
2007, newly reprocessed in 2013, including 2006-2007 data not included in the canonical 
Wittd (120041) analysis. Our Endings are summarized as follows: 


1. We first perform a purely empirical analysis of the He beam observed by Ulysses, using 
2-D Gaussian fits to quantify beam location, amplitude, and width; in order to see how 
these parameters vary with orbital phase and time. There is no evidence for intrinsic 
variability in beam direction by more than a few tenths of a degree. The beam widths 
are also invariant within measurement uncertainties. 
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2. The lack of variability in the He beam from 1994-2007 makes it harder to explain the 

Ulysses discrepancy in terms of He flow variability, with less than a two-year 
time gap for that change in flow to occur before IBEX observations begin in 2009. 

3. We demonstrate that the most powerful constraint that the Ulysses/GAS data provide 

on the He flow vector lies in the shape of the horseshoe-shaped track made by the 
He beam over the course of a Ulysses orbit, which breaks parameter degeneracies that 
plague attempts to derive He flow parameters from individual beam maps. Only by 
htting the Ulysses/GAS data collectively can this diagnostic power be fully utilized. 

4. From a global £t to the full Ulysses/GAS data set from 1994-2007, we hnd the following 

He flow parameters: Vjsm = 26.08±0.21 km s“^, A = 75.54±0.19°, (3 = —5.44±0.24°, 
and T = 7260 ± 270 K. Separate global hts to the three separate fast-latitude scans 
confirm a lack of He flow parameter variability, with the flow direction varying by no 
more than ~ 0.3° and the velocity by no more than ~ 0.3 km s“^. 


5. O ur (VjrmAS/ ) He vector is consistent with the p revious canonical me asurements of 


Wittel (120041) and the contemporaneous analysis of iBzowski et al.l fl2014j ). proving that 


the IBEX/ Ulysses discrepancy is not a product of different analysis approaches. The 
cause of the IBEX/ Ulysses discrepancy remains unknown. 

6 . Our temperature meas urement (T = 7260 ± 270 K) is somewhat higher than the T = 
63 00 ±340 K valu e of Witte ( 2004 ) , but is consistent wit h the T = 7500 K best-£t value 
of iBzowski et al.l (120141) , although iBzowski et al.l (120141) quote a much larger range of 
acceptable values (5500 < T < 9000 K). A higher temperature potentially mitigates 
the IBEX/ Ulysses discrepancy somewhat, as the parameter dependencies in the IBEX 
analysis are such that a higher T would lead to lower A and higher Vjsm, moving the 
IBEX-derived values closer to those of Ulysses. 


7. Using a technique to infer photoionization loss rates solely from Ulysses constraints, we 
compute He densities and loss rates for the three epochs of observation (1994-1996, 
2000-2002, and 2006-2007). Our loss rates agree well with expectations from solar EUV 
constraints for epochs 2 and 3, but not for epoch 1, where we seem to overestimate the 
rate. Our hnal de nsity rn e asure ment is riHe = 0.0196 ± 0.0033 cm“^, somewhat higher 
than estimated by Witte ( 2004 ). riHe = 0.015 ± 0.003 cm“^. 


We would like to thank Dr. E. Mobius for useful discussion and suggestions. This work 
has been supported by NASA award NNH13AV19I to the Naval Research Laboratory. 
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Table 1. 


Ulysses Measurements of the ISM He Flow 


Years 

Fit Type 

VlSM 
(km s“^) 

A (J2000) 
(deg) 

/3 (J2000) 
(deg) 

T 

(K) 

Source 

1994-2007 

Global Fit 

26.08 ±0.21 

75.54 ±0.19 

-5.44 ± 0.24 

7260 ± 270 

this work 

1994-2007 

Single Maps 

(26.08) 

75.60 ± 0.34 

-5.38 ±0.37 

7090 ± 370 

this work 

1994-1996 

Global Fit 

26.14 ±0.24 

75.53 ±0.25 

-5.50 ±0.29 

7230 ± 360 

this work 

2000-2002 

Global Fit 

26.08 ±0.23 

75.57 ±0.22 

-5.39 ±0.27 

7230 ± 310 

this work 

2006-2007 

Global Fit 

25.92 ±0.25 

75.42 ±0.17 

-5.30 ±0.30 

7270 ± 220 

this work 

1990-2002 

Single Maps 

26.3 ±0.4 

75.4 ±0.5 

-5.2 ±0.2 

6300 ± 340 

Witte (2004) 


Table 2. Density and Photoionization 


Years 

UHe from Fit 

(3 ion from Fit 

Aon from EUV'^ 


(10“^ cm“^) 

(10-^ s-b 

o 

1 

-4 

1 

1994-1996 

2.28 ±0.14 

1.36 ±0.27 

0.9 

2000-2002 

1.97 ±0.18 

1.79 ±0.31 

1.6 

2006-2007 

■‘-•’J^-0.31 

n c;7±0.92 
' -0.57 

0.7 


^From Bzowski et ah (2013). 
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240 260 280 300 

Longitude 

Fig. 1.— A WFOV helium beam map made by Ulysses/GAS on 2001 January 24, in ecliptic 
coordinates. White dots indicate actual scan positions used to make the map. In addition 
to the broad He beam, several stars are also visible. 


Latitude 



240 250 260 270 280 240 250 260 270 280 240 250 260 270 280 

Longitude Longitude Longitude 


Fig. 2.— The left panel is a WFOV helium beam map from 2001 July 6, in ecliptic coordi¬ 
nates, the middle panel is a 2-D Gaussian £t to the He peak, and the right panel shows the 
residuals of the fit. The count rates shown are in cts s“^ pixel“^. 
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Fig. 3.— (a) The center of the helium beam as measured by the 2-D Gaussian fits, in 
ecliptic coordinates, with different symbols indicating measurements from the three different 
fast latitude scans. During the course of a fast latitude scan the He beam traces out this 
horseshoe pattern in the sky. The solid line is a polynomial fit to this horseshoe. The numbers 
along the horseshoe indicate orbital phase, where an orbital phase of 0.5 corresponds to the 
ecliptic plane crossing near perihelion, (b) Longitudinal discrepancies from the average 
beam location shown in (a), with the dashed line indicating the ecliptic plane crossing, (c) 
Latitudinal discrepancies from the average beam location, (d) Total angular offset from the 
average beam location. 
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Fig. 4.— (a) Peak count rate amplitude of the helium beam based on the 2-D Gaussian 
hts, as a function of orbital phase, where an orbital phase of 0.5 corresponds to the ecliptic 
plane crossing near perihelion (dashed line). Different symbols indicate measurements from 
the three different fast latitude scans. The count rates shown are in cts s“^ pixeD^. The 
other three panels show (b) background count rates, (c) longitudinal beam widths, and (d) 
latitudinal beam widths. 
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Fig. 5.— (a) Detector geometry function for the WFOV (solid line) and NFOV (dotted line) 
detectors, (b) Detection efficiency as a function of incoming particle energy. 
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Fig. 6.— (a) A WFOV helium beam map from 1995 August 20, in ecliptic coordinates. The 
diamond indicates the expected beam position based on the Witte (2004) vector (assuming 
T = 0 K). Yellow dots are scan positions, and white dots are the measurements within 25° 
of the diamond, which are the data considered in fits, (b) A synthetic image representing 
a fit to the image (FIT ^1) in (a), where the Witte (2004) vector was the initial guess 
for the best-£t parameters, (c) Observed versus predicted count rates for the fit from (b). 
Parameters of this fit are listed, as well as the reduced value {xt = 1-06934). (d) Another 
£t (FIT ^2) to the data with a very different initial guess for the fit parameters from (b). 
(e) Analogous to (c) but for FIT #2. Note the very different best-£t parameters compared 
to (c). 
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Fig. 7.— The solid line in each panel indicates the variation of He beam location with orbital 
phase expected from the Witte (2004) flow vector. Also shown in each panel is what happens 
to this horseshoe-shaped beam track when the latitnde (/9), longitnde (A), or velocity (Vism) 
of the vector are increased (dashed lines) or decreased (dotted lines) by the amonnt indicated. 
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Fig. 8.— (a) Inferred He densities from a global fit to 238 He beam maps, with no correction 
for photoionization losses. The phase dependence is indicative of the presence of such losses. 
Different symbols are for different epochs, as in Figures 3-4. (b) Background count rates 
measured for the 238 maps in the global £t. 
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Fig. 9.— Helium flow parameter measurements from a global fit to the Ulysses /GAS data. In 
each panel, we plot as a function of one of the four He flow parameters of interest, where 
each point represents a separate £t with that parameter held constant and the other three 
parameters allowed to vary freely. Solid lines show polynomial fits to the data points. The 
horizontal dashed line corresponds to a 3a contour, which we use to dehne our uncertainty 
range in each parameter (shaded regions). Vertical dotted lines are the Witte (2004) results, 
which are only signihcantly discrepant for the temperature. 
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Fig. 10.— The top row shows a He beam map from 1994 November 9, a simulated map based 
on our best global £t to the Ulysses data, and the residuals of this fit. The bottom row is 
a similar set of images for a He beam map from 2007 November 26. The diamond indicates 
the expected beam position based on the Witte (2004) vector (assuming T = 0 K). Yellow 
dots are scan positions, and white dots are the measurements within 25° of the diamond, 
which are the data considered in the £t. 
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Fig. 11.— The solid line is the observed track of the He beam center dnring the Ulysses orbit, 
from Figure 3a, which is compared with the expected beam center based on our best-£t He 
flow vector (dotted line) and the IBEX-derived flow vector (dashed line) from McComas et 
al. (2012). The expected beam centers are computed assuming T = 0 K, and the difference 
between the best-£t and observed tracks indicates the level of uncertainty introduced by 
that assumption (see text). Dots mark orbital phases between 0 = 0.30 and 0 = 0.85 in 
increments of A0 = 0.05. 
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Fig. 12.— Helium flow parameter measurements from global fits to Ulysses/GAS data from 
the three fast latitude scans in 1994-1996, 2000-2002, and 2006-2007. The shaded regions 
correspond to the error bars dehned in Figure 9. The horizontal dotted lines are the Witte 
(2004) parameters. We see no evidence for any He flow parameter variability from 1994-2007. 
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Fig. 13.— (a) Interstellar He densities computed from the second fast latitude scan (2000- 
2002), plotted versus orbital phase, assuming photoionization loss rates of = 0 s“^ (plus 
signs), 1.79 X 10“’^ s“^ (circles), and 3.0 x 10“’^ (X’s) s“^. The f^ion = 1-79 x 10“^ s“^ case 
is where the variability of the densities is minimized, (b) Inferred He densities plotted as 
a function of orbital phase for the hrst (diamonds), second (circles), and third (squares) 
fast latitude scans. For each of the three sets of measurements the loss rate is chosen to 
minimize variability in density with orbital phase. These rates are 1.36 x 10“^, 1.79 x 10“^, 
and 5.7 x 10“® s“^ for scans one, two, and three, respectively. 




